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Just as any state of a single qubit or 2-level system can be obtained from any other state by a 
rotation operator parametrized by three real Euler angles, we show how any state of an n-qubit or 
2"-level system can be obtained from any other by a compact unitary transformation with 2"^^ — 1 
real angles, 2" of which are azimuthal-like and the rest polar-like. The results follow from a modeling 
of the Hilbert space of n-qubits by a minimal left ideal of an associative algebra. This representation 
is expected to be useful in the design of new compact control techniques or more efficient algorithms 
in quantum computing. 



PACS numbers: 03.67.-a, 03.65.Fd, 03.67.Mn, 02.20.Sv 

The general control of an A^-level quantum system is 
a subject of considerable interest in many fields, includ- 
ing chemical dynamics and quantum computing. A gen- 
eral statement that establishes the controllability basi- 
cally says that for an A^-level system to be completely 
controllable, it is sufficient that the free-evolution Hamil- 
tonian, along with the interaction Hamiltonian (which 
could involve a sequence of steps) and all possible com- 
mutators among them, form a Lie algebra, which ingen- 
eral is u{N), with iV^ independent real elements 
Once quantum controllability is established, it is impor- 
tant to find the optimal method to actually control the 
system. This involves optimization of both resources and 
time (for examples, see IMQI^ and references therein). 

In this letter, we show that an arbitrary state of an A^- 
level quantum system can be generated from any other 
state by a compact, rotation- like unitary transformation 
with 2N — 1 real parameters. We explicitly demonstrate 
this controllability for n-qubit systems where N — 2". 
We show that a complex algebra of only N indepen- 
dent elements is sufficient to describe an evolution from 
any given reference state to an arbitrary state of the 
system. The N — 2^ independent elements of the lin- 
ear space of the algebra span the iV-dimensional Hilbert 
space, and the A'^ — 1 independent elements excluding 
unity generate polar-like rotations. When combined with 
the A^ independent phase angles for each term, which 
arise from azimuthal-like rotations, there are a total of 
2"+^ — 1 — 2N— 1 real control parameters for the system. 
Thus, it is shown that the minimum number 2A^ — 1 of 
real parameters needed to identify an arbitrary state is 
also sufficient to create it from any other state by a com- 
pact, rotation- like unitary operation. Our results gener- 
alize the well-known description of a single-qubit system, 
any state of which can be specified by the three Euler 
angles of a rotation that relate the state to a "spin-up" 
reference state. 

Clifford's geometric algebra has been used to study Lie 
groups without the use of matrices ^ . The Clifford alge- 
bra C£n of an A^-dimcnsional Euclidean space contains 
a scalar, A^ linearly independent vectors or directions, 
(^) ~ ^ (-^ ~ ^) /■^ linearly independent bivectors or 



planes, (^) linearly independent trivectors, and so on for 
a total of 2^ linearly independent elements. Its unimodu- 
lar even elements, called rotors, form the Spin(A^) group, 
which is a two- fold cover of the rotation group SO(A^) 
and is generated by its bivectors. The 2-state spinor of 
a single qubit can be represented as the projection of 
a rotor in the Clifford algebra Cis of 3-dimensional Eu- 
clidean space onto a minimal left ideal. 0) S Isl • An 
equivalent representation has been extended to the Lie 
algebra of two-qubit systems • This letter shows how 
this structure can be simply implemented for an n-qubit 
system, and more generally, for an A^-level quantum sys- 
tem. 

The qubit or 2-state quantum system is well known 
in many contexts, including the Bloch sphere for spin- 
1/2 magnetic resonance J^] and the Feynman-Vernon- 
Hellwarth model 0| , both of which ultimately have the 
same mathematical structure. The states of a single 
qubit can be represented by elements of Spin(3), which 
is isomorphic to SU(2), projected onto a minimal left 
ideal of €£3 by the projector P3 = i (1 -t- 0-3) with 0-3 a 
unit vector. Projectors are real idempotent elements that 
obey P3 — — (J3P3 ■ In a matrix representation of Ci^, 
the unit vectors aj can be represented by the Pauli-spin 
matrices, cto = 1 is represented by the unit matrix, and 
the volume element is cf 1(72(73 = il. The projector acts 
from the right to reduce the matrix representation of the 
rotor i? to a spinor ip — RP3 with a single non-vanishing 
column. Within the minimal left ideal, the basis states 
of the two-level system can be taken as 

^^3 = fT3P3 ("spin up") (1) 
(T1P3 = cri3P3 ("spin down") 

where (J13 = exp (7rcri3/2) is the bivector cis = 
(Ti(T3 = — CT3cri that generates rotations with rotor R — 
exp {6(713/2) in the plane of ai and (T3, and we noted that 
(7^ = 1 = — 17^3. A general state is a linear combination 

(co + ciai)P3 (2) 

with complex scalar coefficients cq, ci. If the state is nor- 
malized: |co|^ + |ci|^ = 1, we can write cq = e'"cos6'/2 



2 



and ci = e*'" sin 0/2 with real phase angles a,P, and ex- 
press its spinor by 



e cos - 
2 



e'^aissin^ ] P3 



(3) 



= e'" (^cos - - e-'fai3 sin - j P3, = « - /? 

The left factor multiplying P3 is here the rotor for a rota- 
tion with Euler angles {(j>, 9, x)- Its three real parameters 
are the polar angle 6 and two azimuthal angles x ^-nd 4'- 
In terms of physical operators, cj) can be a phase gener- 
ated by free evolution, a phase generated by an applied 
field, and x one generated by both free evolution and 
geometry. 

In this paper, we extend this formalism to systems of 
N states, and represent such states by unitary elements 
of a compact rotor-like group with 2N — 1 real param- 
eters (angles), projected onto a minimal left ideal, and, 
most significant for the control of the system, we show 
that transformations between arbitrary states are given 
simply by such rotor-like operators. 

The proof involves a simple algebraic representation of 
the states. In the case of a single qubit, any state is given 
by a direction in three-dimensional space (on the Bloch 
sphere) together with an overall complex phase. However 
as seen in , it is more simply represented by a complex 
associative algebra Ai of the left factor with just one spa- 
tial dimension times a projector P3 outside the algebra. 
Every element of Ai for a single qubit is a complex linear 
combination of {1, ei} and Ai can be identified with the 
Clifford algebra C£i. We show below that an arbitrary 
state of a system of n qubits can similarly be represented 
by the minimal projection of a left-factor algebra An of 
N = 2^ independent elements generated by products of 
n fundamental ones. The algebra An can be identified 
as the Clifford algebra Cln over the complex field, but 
other algebras are also possible and it may be simpler to 
use an abelian algebra. In all cases, the left factor can 
be written as the product of compact rotor-like unitary 
transformations that can be factored into — 1 polar- 
like rotors and N azimuthal-like ones. The generators of 
the polar-like rotations give the basis states of the sys- 
tem and the polar-like angles themselves give the ampli- 
tudes of the basis states in the state expansion, whereas 
the azimuthal-like angles determine the complex phases 
of the basis-state expansion and are related to the geo- 
metric phases that can be found for A^-level systems by 
an interative procedure recently discussed by Uskov and 
The total number of real angular parameters is 
thus exactly the minimum 2N — 1 parameters needed to 
define any state of the system. 

The Hilbert space of an n-qubit system is spanned by 
tensor products of n single-qubit states. A complete basis 
of states is {|£) = b (^) ^p, £ = 0, . . . , TV - 1} , where the 



reference state ip^ = P = P-i® Pz® ■ ■ ■ ® P-i is &. primitive 
projector, and b {(.) are the hermitian operators 



A,.-i 



A„_2 



'1 ' 



(4) 



is or 1, with corre- 



b {(.) = ® ai 

where each Aj, j = 1,2, 
spending here to "spin up" and 1 to "spin down" . To 
be explicit, we take Xj to be the bits in the binary ex- 
pression of £. In particular, b (0) = 1 because the binary 
number has all bits = 0. Thus, we denote a general 
state ijj by 



N~l 



N-l 



^P^J2 K>Q- ^Qb(£)^o- 



(5) 



«=0 



The b{i) , £ £ [0, n] nZ, form a closed abelian group with 
[b {i)Y = 1. The inner product of states ip and ip is given 
algebraically by {ip'\^) = N (^ipip'^^ , where (x)^ is the 
scalar part of x, that is, the part that is scalar in each 
factor of its tensor product. It is easily verified that the 
basis states \£) are orthonormal. Since ip is assumed to 
be normalized, — J2e=o^ = 1- A direct power- 

series expansion verifies that 



exp (— ib6'/2) = cos - — ?b sin -, 



(6) 



and a product of iV — 1 such compact unitary factors for 
some set of — 1 real "polar angles" Bg gives 



n ''^p 



--b{£)e, 



exp 



N-l 



N-l 



1=0 



with complex coefficients normalized so that 

N-l 



E 

1=0 



1. 



(7) 



(8) 



The expression ((TJ is similar to the state expansion jSJ, 
but the phases are not necessarily the same. To adjust 
the phases, we include N "azimuthal angles" in rotations 
generated by 



ctI\ t£ [0,n]nZ. (9) 



The operators 3 {£) are analogous to the b {£) but with 
CT3 replacing ai. Each 3 {£) either commutes or anticom- 
mutes with any given b {£') and is projected to unity by 
= P. Thus, a product of rotors 



N-l 

R = W exp 



t=o 



exp 



3 (^) 4>t 



2 

N-l 



N-l 



1=0 



t=0 



exp 



b (^') e,, 
Y,b{£')e, 



2 

N-l 



(10) 



3 



when projected onto the minimal left ideal by ■0q gives a 
unique state V' = R'4'o © with arbitrary phase factors. 
Again in terms of physical operators, the 0's are phases 
that arise due to applied fields and 0's arise mainly due 
to free evolution, but sometimes due to geometry. 

The N basis-state operators b {£) are not unique. We 
can alternatively choose elements generating the Clifford 
algebra C£„ with 



ei = 1 • • • C 

e2 = 1 (g) • • • C 
e„, = 6(2"-! 



)3(2 



1 «) (71 = 



b (1)3(0) 
= b (2)3(1) 



(11) 



m—l 



1 



These anticommute and square to +1 : ejek + ej-ej = 
26jk ■ Many other bases are possible in the form 
= b{i)}{Lii) for the expansion (0. Since 
3 (i^) V'o = V'oi tl^s basis states are actually the 
same, but the difference in operators is important since 
given any set {23 (i)}, one can find linear combinations 
X^fco^ qS {£) that are singular, that is noninvertible, and 
such nonunitary cannot be expressed as a rotor. Such 
cases do not cause problems since other operator sets 
can be found in which the linear combinations are not 
singular. 

The strength of our approach is evident from the ex- 
ample of the two-qubit or four-level system rewritten in 
standard matrix form (although the algebra does not re- 
quire matrices) . Four-level systems have been previously 
studied in the context of coupled qubits 
spin 3/2 fermions and controllability |2d|2l|. The 
general state is (jSJ with N — i and 



% = P = P3 



1 



>^3 = ^E3(^) 



b (0) = 3 (0) = 1 ® 1 = 



/ 1 \ 

10 

10 

V 1 / 



/ 1 0\ 





\ / 



b(l) 



b(2) 



1(E) ai 



/O 1 0\ 

10 

1 

V 1 / 
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1\ 
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, b(3) = 








1 





1 














1 










1 





0^ 












0/ 



The matrices for 3 {£) are diagonal and are similarly ob- 
tained. Each b {£) and 3 {£) is its own inverse, but linear 
combinations are not necessarily invertible as we demon- 



strate below. However, the state ip is arbitrary: 

/ Co Ci C2 C3 \ 

Cl Co C3 C2 

C2 C3 Co Cl 

V C3 C2 Cl Co J 



V'o 



In terms of compact unitary factors 



/co \ 
ci 
C2 ■ 

V C3 ; 

(12) 



^p = exp 



exp 



J29,,b{£') 



^0- 



We can relate the angles dg' , (pi to the complex coefficients 
Ci by explicit expansions ^ while noting relations such 
as b(l)b(2) = b(3): 



exp 



. 3 

b' = l 



,b{h') 



(71 (72 (73 

COS cos COS 

2 2 2 



"1 (72 (73 

i sm — sm — sm — 
2 2 2 



. Oi 62 O3 
sm — COS — cos — 
2 2 2 



Ol . (72 . (73 

I COS — sm — sm — 
2 2 2 



b(0) 
b(l) 



with the b (2) and b (3) coefficients obtained from that for 
b (1) by cyclic permutation. This, plus the normalization 
condition (jSJ is sufficient to relate the three polar angles 
6j, j = 1,2,3, to the magnitudes of the complex coeffi- 
cients Ci. The pattern of commuting and anticommuting 
pairs 3 {£) b {£') plus the relation 3 (£) P = P, implies 



exp 



with 



^E^^'3(0 



'=0 



h {£) P 



b{£)exp ( --ai ] P 



ao 



ai = cpQ - cj}^ + cp^ - 03 

"2 = cj)Q + Cj)-^ - Cj)2 - (j)^ 

and this is enough to relate the complex phases of the Ci 
to the four azimuthal angles 4>j^, fc = 0, 1, 2, 3. 

As mentioned above, for any choice of basis operators 
*8 {£) , there exist some states for which the linear combi- 
nation c^S {£) is singular, that is, has no inverse. For 
the basis operators b {£) , for example, the problem arises 
with the entangled states V'i = ('' (0) i b (3)) V'o -H 

is seen by inspection that the matrix Y^^=o'^^^ (^) given 
in (|12() has a vanishing determinant when C3 = ±co and 
ci = = C2. The problem is resolved by using differ- 
ent basis operators. Thus, if we adopt the €£2 basis with 
phases that make all basis elements antihermitian, we get 
the expansion 

/ Co ici ic2 -C3 ^ 
ici Co C3 -ic2 
ic2 -C3 Co ici 

\ C3 -iC2 ici Co J 



3 

E 



qQ5 {£) = 



(13) 



4 



which is unitary for the choice cq = C3 = l/\/2 and 

Ci = = C2. 

The method can be appHed to arbitrary 2"-level sys- 
tems. The relation of the tensor-product basis states, in- 
terpreted as giving spin orientations of separate spin-1/2 
systems, to the eigenstates of the system depends on the 
Hamiltonian. For example, for 71 = 2, the Hamiltonian 

H = Eo+ 7iSi • B «) 1 + 1 (g) 72S2 • B -I- 2Asi • 83, (14) 

where B is the external magnetic field, 7^ are the gy- 
romagnetic ratios, and 2A is the spin-spin coupling con- 
stant, has the matrix form 






V 





X 




hujn 



2 / 



(15) 

where a;o±± = ujq ± uJi ± lo2- States ipQ — P and 
tp^ — b (3) '00 = CTs a^P are eigenstates, but linear 
combinations oi ipi — b (1) P and ip2 = b (2) P form the 
other eigenstates. In the 2-D basis {^Aij "02} > 



H = hujo- A/2 + 



A 
A 



A 

-A 



(16) 



with A = h{uji ~ LU2) and the eigenenergies are hioa — 
A/2±\/A2' + A .Since we have a total of four independent 
energy parameters, we expect to be able to represent an 
arbitrary 4-level system with these coupled qubits. 

More generally, however, there is no need to associate 
the tensor-product states literally with independent spins 
on spin-i particles. Tensor products of Pauli spin ma- 
trices are simply a way to specify matrix representations 
of orthogonal states and operators on them. A complete 
set of 2" X 2" matrices is always given by tensor products 
of n spin matrices afj,. Nothing inherent in our treatment 
restricts it to systems of n-qubits. Any set of N levels 
can be similarly handled. In the same vein, the refer- 
ence state tpQ need not be one with maximum total spin 
component along the z axis. Formally, we can trans- 
form the reference state multiplying our spinors ■)/' by a 
transformation rotor from the right. This shows that the 
transformation between any two states of the system can 
be represented by a unitary rotor with — 1 polar-like 
angles and A^ azimuthal ones. For A^ > 1, this is fewer 
than the A^^ real parameters required to specify a uni- 
tary transformation on the system in general, but this 
reduction simply reflects the fact that only a small sub- 
set of the possible unitary transformations are required 
to transform between arbitrary states. For most prac- 
tical problems in quantum control, it is only necessary 
to transform a few given states to others (for example, 
making a quantum gate), and therefore our representa- 
tion can be expected to lead to more compact control 
schemes. 



In summary, we have demonstrated that evolution of 
an arbitrary n-qubit or N = 2"-level quantum state into 
another state can be described by a complex algebra of 
only A^ independent elements. This compact unitary 
transformation consists of 2"+^ — 1 real angles, 2" of 
which are azimuthal-like and 2" — 1 that are polar-like. 
These results are expected to lead to optimal schemes for 
coherent control and quantum computing. 
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